arXiv:chao-dyn/9706024vl 24 Jun 1997 


Quantitative study of scars in the boundary section of the stadium billiard 

Fernando P. Simonotti, Eduardo Vergini, Marcos Saraceno 
Departamento de Fisica, Comision Nacional de Energia Atomica, Av. Libertador 8250, 1429 Buenos Aires, Argentina. 

We construct a semiclassically invariant function on the boundary of the billiard, taken as the 
Poincare section in Birkhoff coordinates, based on periodic orbit information, as an ansatz for the 
normal derivative of the eigenfunction. Defining an appropriate scalar product on the section, we 
can compute the scar intensity of a given periodic orbit on an eigensate, as the overlap beetween 
the constructed function and the normal derivative on the section of the eigenstate. In this way, we 
are able to investigate how periodic orbits scar the spectrum and how a given eigenstate decompose 
into scar functions. We use this scheme on the Bunimovich stadium. 

PACS Numbers: 03.65.Sq, 05.45.+b, 03.40.Kf 


I. INTRODUCTION 

Since the observation of imprints of periodic orbits in quantum eigenfunctions, scars, by McDonald and Kaufmann 
[Ql; a vast amount of work has been done towards the understanding of this phenomenon. The numerical work and 
theoretical analysis of Heller has been of great importance. Bogomolny pushed the theory of scars further, his 
developments relied on the smearing of the probability density over a small energy range. A similar approach, but in 
phase space rather than in coordinate space, was used by Berry Q. A theory for individual eigenstates was developed 
by Agam and Fishman who constructed a semiclassical Wigner function. The integration of this Wigner function 
in a narrow tube along a periodic orbit gave them the scar weigth. Smilansky Q used the scattering approach to 
define a function on the Poincare section, which was tested for scars. An important tool in the following analysis, the 
stellar representation, was developed by Tualle and Voros 0. 

In this paper, we construct a semiclassicaly invariant function on the Poincare section, built on a given periodic 
orbit, which we call the scar function, which can be extended to the domain via the Green theorem (D- We define 
the scar intensity as the overlap beetween this scar function and the corresponding reduction of the eigenfunction on 
the section, with a given measure, so as to mimic the overlap in the domain. Using this construction on the stadium 
billiard, we are able, by means of the symbolic dynamics, to identify scars of single periodic orbits and of families of 
them in the quantum spectrum. Also, as we do not resort to energy smearing, we can decompose an eigenstate in 
periodic orbit functions (the scar functions.) 

This paper is organized as follows. In section || we formulate our approach to the calculation of a scar measure, and 
introduce the necessary objects. In section g we compute scars for the Bunimovich stadium and analize how different 
families of periodic orbits scar the corresponding eigenfunctions, particularly the whispering gallery and bouncing ball 
families. Our conclusions and closing remarks can be found in the last section. 


II. THE SCAR FUNCTION ON THE BOUNDARY 

All the information of a given eigenfunction of the billiard (with Dirichlet boundary conditions) is contained in 
its normal derivative evaluated on the boundary. By means of the Green function we can obtain the wave function 
'l'y(r), with wave number in the domain in terms of its normal derivative: 

^'^(r) = ds Go{K-, r, r'(s)) ^[^(r'(s)) (1) 

where Go(fci,;r,r'(s)) = {—i/4)HQ^\k,,\r — r'(s)|) is the free Green function |^. is the Hankel function of the 

first kind. 
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We thus treat the normal derivative as the fundamental irreducible object that will be tested for the presence of 
scars. The function 


=-^(r'(s)) ( 2 ) 

can be thought of as the coordinate representation of an abstract vector in the Hilbert space of periodic square 
integrable functions on the boundary. This boundary is also the standard Poincare section for the classical dynamics 
and reduces the motion to a canonical mapping in the Birkhoff coordinates (g,p) 0. The coordinate g is related to 
the arc length coordinate at the boundary where the bounce takes place by g = (s/perimeterj^^p^j and p = p-t/|p| 
is the fraction of tangential momentum at this point. 

The Fourier transform of (|^) would represent it in the momentum representation and a coherent state one would 
display its features in the Birkhoff coordinate plane. We follo w to choose this route. Other representations are, of 
course, possible and have been used in this context. Smilansky H, for example, has used angle and angular momentum 
variables as phase space coordinates, which is the natural description when employing scattering methods. However, 
in this basis the representation of diffraction effects on functions on the boundary is singular and difficult to observe. 
The Birkhoff coordinates, besides embodying the natural geometry of the billiard, avoid this problem by prescribing 
definite periodic functions as candidates for boundary eigenfunctions. 

As shown in |lH , for eigenfunctions normalised to unity in the domain, in a neighbourhood of a given k value, we 
have the quasi-orthogonality relation for the normal derivatives 

ds n ■ r 0(1). (3) 

Thus, with this measure, the set of eigenfunction in a narrow range of k is orthonormal and span a linear space of 
dimension 0{k). 

It is then convenient, as we want to work exclusively on the boundary, to adopt a definition of scalar product 

((<('1^)) = ^ / dsn-v (4) 


Any of these functions can be extended to the domain by means of Green’s theorem; using it as an ansatz for the 
normal derivative, ff-(s), and setting 4'(s) to zero. 





F)C^ 


Go(k;r,r'(s)) |^(r'(s)) 


(5) 


Of course, is not an eigensolution, because the limiting value of d>(r), as r goes to the boundary, is not zero; that 
is to say, the function is discontinuous at the boundary. This function depends on the continuous parameter k, which 
controls the semiclassical limit. 

With this definition, the normal derivatives 4’v(s) of eigenfunctions normalized in the domain are orthonormal (to 
order 1/k) in a small range Ak = (2 x perimeter/area). 

For the phase space representation we construct coherent states with the correct space periodicity |Q , defined as 


3\p g) = 


-) 

aiT J 


1/4 oo 

exp {i k p {s — a)) exp 




( 6 ) 


This is a boundary wave packet, periodic in s, which is localized at the point (p, g) in the Birkhoff Poincare section 
phase space. 

A single wave packet represents a bounce off a specified point on the boundary with a given tangential momentum. 
Thus, to extract the phase space contents of a given eigenfunction, we can construct the overlap 
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Thus, a first visual display of the eventual localization and scarring of the eigenstates comes through the Husimi 
function 


T-i-Ap.q) = \Au{p,q)\ ■ 

We show an example of this for the stadium eigenfunction with k = 100.954920427642 in Figure ^ 
Clearly, a first quantitative measure of the scarring of a periodic orbit would be 

1 ^ 

l) = (Pli - 97 .) 

2=1 


( 8 ) 


(9) 


which averages the probability over the N points 7 ^ of the periodic trajectory 7 . This measure was used by Muller 
and Wintgen in the context of the diamagnetic Kepler problem |^. 

This average over the probabilities of the periodic points does not take into account the phase relations, due to 
semiclassical propagation, beetween them. Therefore, it seems a more convenient strategy to average the amplitudes 
with the proper phase differences. A better measure, then, is provided by the construction of the scar function 


(s|(/ 7 (fc, 7 )) 


1 

K(</5(fc,7)|(^(fc,7)))|V2 


N 

{s\pjqj} 

i=i 


( 10 ) 


where qi and pi are the Birkhoff coordinates of the periodic points. 
The phases fj are defined by 


fj — klj jjT ^ 17 


( 11 ) 


where Ij is the distance in configuration space beetween the initial point of the periodic orbit and the j — th point. 
The second term takes into account the boundary conditions (Dirichlet) and the third, the conjugate points along the 
trajectory. The inclusion of these phases is very important in the determination of the existence of scars. 

The total accumulated phase, /jv, will not be, in general, a multiple of 27r. In order to have an invariant function, 
depending only on the orbit and not on the starting point, we add an additional phase to each point, so as to make 
/jv = 27rn, with n an integer: 


fj - fj + ( 12 ) 

with a the minimum beetween (/Ar)jjiod 2 tt 27r — (/Af)niod 27 r' 

This state is a coherent sum over a periodic orbit and, thus, is a good candidate for an invariant probe depending 
only on the orbit. 


A. The Scar Intensity S^{kA and the Scar Length Spectrum sAi) 

We define the scar intensity 

S{v,l) = (13) 

Notice that k is set to ki, in {p{k,^)\. This measure of the scar intensity differs from mainly by interference terms. 

Each wave packet in ( 0 ) represents a localized plane wave hitting the boundary at a specified point in a specihed 
direction. Thus (0), when seen in this light, can be assimilated to a superposition of plane waves which privileges 
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the wave directions associated to the periodic orbit. For example, in Figure ^ we show a periodic orbit of the stadium 
billiard and the associated scar wave function in different representations, with k = 100.954920427642 (the same as 
in Figure |^). 

The scar wave function in the domain is a solution of the Helmholtz equation with a given value of k. So, if 
we expand it in terms of the exact eigenfunctions, we expect that the more significant contributions come from the 
eigenfunctions with closer k^ to k (more precisely, \k — k^\ < (2 x perimeter/area). Then, using the quasi-orthogonality 
relation i) the norm of the scar wave function in the domain is 1 to order k ^. Our aim is to describe this subspace 
in terms of states constructed on periodic orbits, as in (p^. This is not dissimilar to the task of describing them in 
terms of plane or cylindrical waves. However, the peculiar linear combinations taken in (10), being semiclassically 
invariant under the bounce map, should provide the most important correlations. 

The density of states of the billiard has a semiclassical representation as 


dsciik) « (d(fc)} 


OO 

E E 

p r—l 


|det(I-Tf)|(i/2) 


cos {r{klp — r'p 7 r/ 2 )) 


(14) 


where the first sum is done over primitive periodic orbits, the second sum takes into account their repetitions; Ip is 
the length of the orbit, i/p is the Maslov index, Tp is the monodromy matrix. The smooth part of the density of states 
is given by {d{k)). 

The Fourier transform of 0 ) provides a distribution linked more directely to the classical motion, i. e., the length 
spectrum, which shows well defined peaks at the lengths (actions) of periodic orbits. 

In order to focus more specifically on the scarring features of a single orbit along the spectrum, we set 7 to a given 
periodic orbit in (|l^ ) 

S^{k,,) = \{{ipj{k^)\(l},,))\‘^. (15) 

Its Fourier transform is 

s,(i)=y: S-si{kv) exp (ifcj/Z), (16) 

ku 


which we will call the scar length spectrum. 

If periodic orbits obeyed Bohr-Sommerfeld like quantization conditions, we would expect a periodic behaviour of 
S-y{k,^) with a period Afc = This implies periodic sequences of scarred states along the spectrum . These 

states have been also observed in other billiards [| §• Here we test for these periodicities directly in the scar length 
spectrum. The periodicities would be exact but for the fact that the periodic orbit basis is not orthonormal, and the 
quantization rule of a single orbit does not lead necessarily to a quantized state. 

To eliminate spurious behaviour of the Fourier transform due to end effects of the /c-interval, we multiply the scar 
intensity in 0 ) by a function vanishing at the ends; typically a quadratic function p^ . The scar length spectrum 
still shows large fluctuations with very clear local average peaks at certain lengths. This local structure is evidenced 
by averaging the resulting Fourier transform. 


III. THE BUNIMOVICH STADIUM 

The boundary of the stadium billiard is defined by two semicircunferences connected by two straight segments. Of 
all the possible stadia, we consider only the one with relation 2 : 1 beetween its total length and height, and we will 
scale the lengths in such a way that the perimeter is 4 -|- 27r and the area is 4 -|- tt. 

At the classical level, we will describe the periodic orbits following Biham and Kvale ||^. Their symbolic dynamics 
is a six-symbol one where each symbol corresponds to a bounce off the boundary: 

• 0: A bounce off the lower straight segment. 
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• 1: A clockwise bounce off the left semicircle or a single anticlockwise bounce off the left semicircle. 

• 2: A bounce off the upper straight segment. 

• 3: A anticlockwise bounce off the right semicircle or a single clockwise bounce off the left semicircle. 

• 4: A not single anticlockwise bounce off the left semicircle. 

• 5: A not single clockwise bounce off the right semicircle. 

A bounce is a single one if it is not preceded or followed by a bounce off the same section of the boundary. 

This dynamics has to be pruned. This pruning is geometrical and corresponds to the symbolic dynamics of the 
stadium of inhnte length. As the length is made hnite, more pruning rules appear (of a dynamical character) as 
described in [l^ . 

Using this symbolic description we have computed all the periodic orbits up to 10 bounces and a few selected ones 
of much higher periods, and we have ordered them (somewhat arbitrarily) by their number of bounces and symbolic 
codes. In Figure || we show the correspondence of these calculated orbits as projections in configuration space and in 
the Birkhoff Poincare (desymmetrized) section. 

At the quantum level, we compute the energy levels and eigenfunctions by the scaling method ||ll|, which gives 
directly all eigenvalues and eigenfunctions very precisely and efficientely. We have computed 1654 consecutive levels 
and their eigenfunctions ranging from k ~ 62.8 and k w 125.2; and other selected ones. 


A. Periodic Orbit Decomposition of Eigenfunctions 

In Figure || we demostrate the advantage of using the present scar function as opposed to the simple unphased 
average. S', of We take one of the scarred wavefunctions of Heller (Figure ^ and plot the quantities S and 
S' as a function of the periodic orbit label 7 (Figure ||). The periodic orbits are ordered in increasing periods and, 
within each period, by symbolic codes. Recurrences in S sometimes occur due to the existence of orbits of period 
n X p that almost retrace n times the periodic orbit of period p. However, many recurrences are also due to short 
homoclinic and heteroclinic excursions. 

The peaks are more clearly defined in the S plot, due to the enhancement brought about by the semiclassical 
dynamics fj (note that the maximum scar value is around 0.114 for S and 0.06 for S'). Due to the binning of the 
interval, the strongest peak of S' overlaps with one of the secondary peaks of S. In the leftmost inset in Figure ^ we 
show the periodic orbit that scars this eigenfunction the most. Moreover, the two most prominent secondary peaks 
come from two different homoclinic excursions of this orbit (see rightmost inset in Figure 

Both measures indicate the presence of scars, i.e., amplitudes larger than the average fluctuation. However, as 
S captures the phase relations of periodic orbits, the basis \ip{k,^)) is ‘closer’ to reflect the invariant properties 
characteristic of the stadium. Thus, we expect the amplitudes S to have much larger fluctuations (and, therefore, 
clearer scars) than S'. This fact shows clearly in Figure]^, where we show the distribution of scar intensities, N{S). 
The small amplitudes are distributed approximately as an exponential. The distribution corresponding to S is much 
broader than that of S'. For this strongly scarred state, there is a large region beetween 0.07 and 0.114 where no 
scar intensities appear. So, \ogiQ{N(S)) goes to —00 in this region. At 0.114 a single periodic orbit gives a large 
scar, yielding \og^Q{N{S)) = 0. The occurrence of this peak is, however, a rare event. The secondary peaks, due 
to homoclinic excursions of the periodic orbit which gives the strongest peak, are to be found at, approximately, 
S = 0.071 and S = 0.066. 

We have tested all 1654 eigenfunctions in the range beetween k Ki 62.8 and k « 125.2 against scarring by the first 
617 periodic orbits (i.e., up to 9 bounces), giving a total of approximately 1 million scar intensities. The distribution 
of these intensities is shown in Figure We can observe three different sections: strong scars, weak scars and the 
region closer to S' = 0, where most of the scar intensities are (there are around a million scar intensities from S = 0 to 
S = 0.1, whereas only 760 are to be found with S > 0.1.) This distribution is very different from the Porter-Thomas 
result. The reason is that the basis of {s\(p{k,'y)) is not orthonormal and, moreover, is chosen so as to be closely 
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related to the dynamics. Should we examine the scars in any other basis, unrelated to the dynamics of the stadium, 
for example a plane wave basis, we would only expect a statistical distribution of the intensities, in accordance to 
random matrix theories. 

The presence of the peaks in the strong scar region in Figure quantifies the scar phenomenon and shows, in 
accordance to Shnirelman’s theorem |p^ , that scarring is exceptional. However, it is the only remaining signature of 
the specific classical behaviour of the system, as embodied in its periodic orbits. 


B. Families of Periodic Orbits 

1. Whispering Gallery Family 

This family is composed by trajectories of code 5"0“1™2^; where a and b are 0 or 1, and m and n are positive 
integers (as there is time reversal symmetry, the same orbit can be described by the code 3"2^4™0“). The whispering 
gallery limit is approached as TO,n —> oo simultaneously. For the periodic orbit to exist, as this limit is approached, 
the difference beetween m and n should remain finite; the largeness of this difference being determined by the length 
of the stadium. This is a clear example of dynamical pruning (as opposed to geometrical pruning, i.e. independent of 
the length of the stadium) 

We have found that the whispering gallery trajectories that show a more pronounced periodicity in the scarring 
are those defined by m > 2 and n > 2, independent of their symmetry or value of a and b. Of course, where high 
symmetry is present, the periodicity is stronger. 

We show in Figure ^ the function Sj{k) and the scar length spectrum for some whispering gallery periodic orbits. 
We see how the scar length spectrum shows clearly the periodicity of S^{k), defined by the fundamental length of the 
orbit and its repetitions. (The lengths shown are multiples of T-y/2 because of the symmetry of the orbit.) 

The width of the group of states that participate in the scarring is constant in the region considered. This means 
that more and more states are involved in one ‘Bohr-Sommerfeld’ interval. However, only a few, typicaly one or two, 
show visible scars. In Figure || (first panel), where the width of the groups is approximately Ak « 1.27, the first group 
involves around 35 states at fc « 62.88 and the last involves 45 at fc w 125.60 in accord to the change in the density 
of states. 

We can exemplify how well the scar intensity picks up the scarred eigenstates. We look for scars in the region 
delimited by fc « 77.14 and k « 77.34 for the periodic orbit with code 55551111 (see Figure P). We find two 
contiguous states that are scarred by this periodic orbit, namely k = 77.22171991174 and k = 77.24033835210. The 
probability densities of these eigenfunctions confirm this fact: see Figure H- Notice the similarity beetween both 
eigenfunctions and beetween them and the mentioned whispering gallery periodic orbit. 


2. Bouncing Ball Family 

This family is composed of three subfamilies, following ||^, whose symbolic codes are 33(02)"11(02)"' (family A), 
23(20)"'21(20)" (family H), 3(20)"1(02)" (family C). The bouncing ball limit is n ^ oo, where the resulting periodic 
orbit has increasingly smaller x component of the wavenumber k. We will consider the first five members of each 
family, with periods ranging from 6 to 24 (see Figure 0 )- 

The bouncing ball eigenstates are approximately described by those of a rectangle with the same size as the one 
inscribed in the stadium. So, the quantized wavenumbers are given by fc « TTy/n^ + where nx{y) is the number of 
nodes along the x(%j) axis. 

The preceding considerations tell us not to expect a scar length spectrum that is peaked in the length of the given 
periodic orbit and its multiples. This is so because there is no such periodicity in the scar intensities for the bouncing 
ball families. This is what we observe for those orbits tending to the bouncing ball limit, see an example in Figure 
[l^. Notice, though, how for big L the peaks appear for even L. This fact is related to the multiple bounces beetween 
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the two straight segments of the billiard, approximately of length 2. Some of the first few orbits of each family show 
a single peak in the length of the orbit, with no peaks (or small ones) in the multiples, see an example in Figure |^. 


3. Other Periodic Orbits 

As ‘scars are scarce’ , the scar length spectrum for most periodic orbits shows no peaks in the associated length. 
Other orbits show rather more complex patterns with many lengths that are not easily assigned to other periodic 
orbits (whether in the homoclinic family of the first or not). However, there is a small set of orbits for which the scar 
intensities have the expected periodicity (some examples in Figure |l4|). 

IV. CONCLUSIONS 

We have constructed a quantitative measure for the presence of scars that, by taking into account semiclassical 
phase correlations, provides a sharper indication of their presence. The measure is constructed as an ansatz for the 
normal derivative of a state representing a pure stationary scar and therefore testing the most the irreducible contents 
of the eigenfunctions as contained in the normal derivative. Thus, it is very suitable for numerical calculations, as 
it involves only boundary integrations, avoding completely- except for graphical display- any integration over the 
domain of the billiard. 

We have provided examples of the decomposition of a single eigenstate into scar functions and of the systematic 
way in which some orbits appear in the k spectrum. The scar intensities show how an eigenstate is distributed on 
a basis of quantum states constructed from periodic orbits. Such a basis is clearly not orthogonal (and, probably, 
overcomplete). The investigation of the properties of this basis remains to be done. 

By testing many orbits and many eigenstates we have found that scars are quite rare, in accordance to expectations 
from Shnirelman’s theorem | ^ . Even less frequent is to find sequences of states scarred periodically (in k) by a given 
orbit. Most eigenfunctions decompose in periodic orbits in such a way that no one prevails over the others. This, in 
turn, implies no clear scars in most instances. 

The families that show stronger scars are the whispering gallery and the bouncing ball ones, both being rather 
exceptional families. 
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FIGURE CAPTIONS 


FIG. 1. Eigenfunction (left panel) and Husimi representation (right panel) for k = 100.954920427642. This function is 
scarred by the periodic orbit with code 23202120. The crosses are the periodic points of this orbit. 


FIG. 2. Different representations of the scar function of a periodic orbit. Upper left panel: Periodic orbit in the fundamental 
domain (symbolic code 23202120, following |^) Upper right panel: Real part of the scar wave function in the boundary. Lower 
left panel: Probability density in the domain (via Green function). Lower right panel: Husimi representation. 


FIG. 3. Periodic orbits of 4 bounces in the configuration space (upper panel) and in phase space (lower panel) In the latter 
plot, the points are depicted by dots or crosses, depending on the sign of p. The dotted vertical line is placed at the value of q 
where the discotinuity in curvature occurs. 


FIG. 4. Odd-odd stadium eigenfunction with k = 130.4886755073. 


FIG. 5. S and S' functions for k=130.4886755073 (upper and lower panel, respectively.) In the leftmost inset we see the 
periodic orbit with maximum scar intensity; in the rightmost one, the orbits that give the two secondary peaks. The orbits are 
ordered by increasing number of bounces, up to 9. 


FIG. 6. Distribution of scar intensities, N{S), for S (full line) and S' (dashed line) for k=130.4886755073. 

The base of the logarithm is 10. 


FIG. 7. Distribution of scar intensities, N{S), for the 1654 consecutive eigenfunctions and first 617 periodic orbits. 

The base of the logarithm is 10. 


FIG. 8. Scar intensity and scar length spectrum for some whispering gallery periodic orbits. Upper left panel: Scar intensity 
for periodic orbit 3332440. Lower left panel: Scar length spectrum for periodic orbit 3332440. Upper right panel: Scar intensity 
for periodic orbit 55511111. Lower right panel: Scar length spectrum for periodic orbit 55511111. 


FIG. 9. Scar intensity for a range of k, for the periodic orbit depicted in the inset. 

FIG. 10. Group of eigenfunctions in the same range of k of Figure 


FIG. 11. Bouncing ball periodic orbits, first five members of familes A, B and G (first, second and third row, repectively) 


FIG. 12. Scar intensity and scar length spectrum for a periodic orbit in the bouncing ball limit. 


FIG. 13. Scar intensity and scar length spectrum for a periodic orbit not in the bouncing ball limit. 
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FIG. 14. Scar intensity and scar length spectrum for some low period periodic orbits. Left panels: Scar intensity and scar 
length spectrum for periodic orbit 3210. Right panels: Scar intensity and scar length spectrum for periodic orbit 2321. 
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